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Abstract 

We show how the phase diagram for spin torque switching in the case of perpendicular anisotropy 
is altered at nonzero temperature. The hysteresis region in which the parallel and antiparallel 
states coexist shrinks, and a new region of telegraph noise appears. In a small sample, the region 
of coexistence of a precessional and parallel state can disappear entirely. We show that the phase 
diagram for both zero and nonzero temperature can be understood and calculated by plotting 
an effective energy as a function of angle. A combinatorial analysis is useful for systematically 
describing the phase diagram. 



I. INTRODUCTION 



In this paper we consider a thin film magnetic element with perpendicular anisotropy. 
The phase diagram for this system has been studied theoretically^ at zero temperature on 
the assumption of a homogeneous single domain, and experimentally 2 -. Some discrepancies 
appear to be due to inhomogeneities, but some may be due to the fact that the zero- 
temperature theory does not take into account thermal fluctuations. The latter effects 
can be calculated semi-analytically within the single domain model, whereas inhomogeneity 
effects probably can only be dealt with numerically. Thus, in this paper we will generalize 
the single-domain phase diagram to nonzero temperature; to the best of our knowledge, this 
has not been done previously. 

In uniaxial symmetry, the energy depends only on the angle 9 of the magnetization from 
the easy axis: 

E(9) = M s [H<f sin 2 9 - H e cos 0] (1) 

where M s is the saturation magnetization, Hffi = — M s is the effective anisotropy, and 
H e is an external field along the easy axis (normal to the film). Any discussion of spin 
torque dynamics begins with the Landau-Lifshitz (LL) equation^ for the time derivative of 
the magnetization, M, sometimes loosely referred to as "torque". The LL equation has a 
precession term which contributes only to the azimuthal component and doesn't change 
the energy; the changes in energy are controlled by the 9 component 

M e = M s 9 = -jsmO{aM 8 [Hf cos 9 + H e \ + J *f s } (2) 

(1 + B cos 9) 

where 7 is the gyromagnetic ratio. The first term (proportional to the LL damping parameter 
a) controls damping and pushes M toward the easy axis, and the second (spin torque) term 
scales with a parameter J proportional to the current (B is Slonczewski's dimensionless 
torque asymmetry parameter-) 

II. EFFECTIVE ENERGY 

Defining an effective energy in the presence of spin torque is nontrivial, since only the 
precession term in the Landau-Lifshitz equation conserves energy; the damping and spin 
torque terms are non-conservative. 
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The most rigorous way to derive an effective energy is by finding a steady state solution 
of the Fokker-Planck equation 3 ' 5 , which has the form exp(— E^/hsT) where 

E eS = M^Hfp (1 - cos 2 9) - H e cos 9 - ln(l + B cos 9)] (3) 

It is worth noting, however, that Eq. [3] can be obtained heuristically; in this special case 
of uniaxial symmetry, if we compute the "work" done by the "torque" (Eq. [2]), the result is 
proportional to the effective energy: J Mgd9 = ajE e Q. 

rr 

We will use a non-dimensional form e e g = -E e ff/ i^oH^M s of the effective energy, 

e eff = \{l-u 2 )-h e u-^\n{l + Bu) (4) 
where u = cos9, h = H e /H^, and j = J/aHj^. 



III. ZERO-TEMPERATURE PHASE DIAGRAM 

The behavior of Eq. H] for various values of field and current is shown in Fig. [TJ The 
only free parameter is the Slonczewski parameter B, which we have taken to be 0.5 in the 
figures. The horizontal and vertical axes are the dimensionless magnetic field and current, 
h and j, and at each point of a 3 x 3 grid there is an inset graph showing the effective 
energy as a function of u = cos 9, at the field and current corresponding to the center of 
the inset. The center graph, at zero field and current, shows only the anisotropy energy 
—u 2 /2, a parabola. At this point there are two stable states (minima of e(u) at u = ±1 - 
note that e(u) need not be flat at minima that lie at the boundaries of the physical region 

— 1 < u = cos 9 < 1). Thus in this region of the phase diagram both the parallel (cos 9 = 1) 
and antiparallel (cos# = —1) states are stable, so we have labeled it "AP & P". 

If we move to the right from the center of Fig. we add the Zeeman term — hu, which 
simply shifts the parabola to the left so its maximum is outside the physical region, and 
there is only one (parallel) minimum at u — 1, and we have crossed a phase boundary (the 
red circles - colors online) to the "P" region where only the parallel state is stable. At 
this boundary, the maximum is just leaving the physical region at u = —1, i.e., de(u = 

— l)/du = 0. On the left, at h = —3.0, the parabola shifts in the opposite direction (right) - 
the maximum passes out of the physical region at the black (square symbols) phase boundary, 
where de(u = l)/du = 0, and only the AP state is stable. 
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FIG. 1: The phase diagram in the h-j plane for the perpendicular spin-torque element described in 
the text, for spin-torque parameter B = 0.5. Lines with symbols are phase boundaries, described in 
text. Insets show e vs. u = cos at a 3 x 3 grid of points at h = —3.0, 0, +3.0 and j = —2.5, 0, +2.5, 
plus one additional inset for the small PS & P region. 

Moving up from the center of Fig. [1] (increasing the scaled current to j = 2.5), the 
effective-energy inset graph includes the logarithmic spin-torque term — ln(l + Bu). This 
has a divergence at u — — 1/B — —2, outside the physical region, but starts to rise at the 
left as seen in the top center inset. This mimics the effect of a positive field (right inset) 
and brings us across the phase boundary to the upper-right single-phase P region. 

At the upper left inset (negative field h = —3.0), the effect of the field (lowering e at the 
left) and the spin torque (raising e at the left) oppose each other, and the spin torque can 
raise the energy at the left enough to create a minimum away from the u — ±1 boundaries, 
physically corresponding to a precessional state - thus this region of the phase diagram 
is labeled "PS". The effective energy in the remaining small sliver of the phase diagram, 
between the curved line (blue triangles) and the straight phase boundaries, is shown as a 
tenth inset - both the precessional and parallel states are stable, so this region is labeled 
"PS & P". 

Making the current negative (bottom center inset) causes a decrease in energy at the 
left side of the inset, mimicking a negative field. The only minimum is now the AP (left) 
minimum, and the system is below the black (squares) phase boundary in the AP region. 
At the lower left (adding a real negative field) further stabilizes the AP state. Moving to 
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the lower right (positive field) restores the u = 1 (parallel) minimum and returns us above 
the black phase boundary to the AP & P region 

IV. SYSTEMATIC COMBINATORIAL CONSTRUCTION OF PHASE DIA- 
GRAM 

It is not obvious from the above approach to the phase diagram how many regions there 

can be. We now enumerate the regions more systematically. The physical behavior depends 

on the number and position of extrema of the effective energy e(u). The condition de/du = 

gives a quadratic equation for u, with two solutions u~ < u + if j < 4(/t — -^) 2 , i.e. below 

the parabola shown by blue triangles in Fig. [TJ Above this parabola there are no flat 

maxima or minima (as opposed to the boundary extrema at u — ±1), and the boundary 

minimum always occurs at u=l: this region is a parallel (P) state. The lowest straight line 

[black squares, j = —{h + 1)/(1 + B)] in Fig. Q] is where de{u = l)/du = 0, i.e. one of 

the extrema occurs at the right boundary, and the other straight line [circles (red online) 

j = — {h— 1)/(1 — B)] has deiu = —l)/du = 0. Below the parabola we classify the positions 

of the extrema relative to the boundaries u = ±1 of the physical region by a four-character 

string such as | x x|, where the vertical lines represent u = —1 and u = +1 respectively 

and the x's represent the minimum and the maximum (in that order). There are exactly 
4 

( ) = 21^7 = 6 ways to pick 2 of the 4 positions for x , i.e. to order these 4 characters; if we 

add the string 1 1 to represent the P region without extrema, these are exactly the 7 light grey 
(yellow online) circles in Fig. El a distorted cartoon of the phase diagram. The boundary 
between each adjoining pair of regions (light circles) is labeled by a darker grey circle (blue 
online) with a string in which two of the symbols are superposed (e.g., * x |, where * indicates 
a | and x superposed, near the bottom center of Fig. El means there is an extremum at the 
left boundary u — — 1; there the | and x pass each other, converting | x x| into x| x |). 
Each boundary similarly involves the crossing of two symbols, except that | and | cannot 
cross, and x and x "cross" only when they merge - beyond this boundary (the distorted 
parabola) there is no minimum or maximum. At the intersections of boundaries (darkest 
grey circles, green online) there are two coincident symbols: at the rightmost intersection, 
x crosses * - at that point the minimum and maximum annihilate at u — — 1. 



5 



FIG. 2: Cartoon version of the phase diagram (Fig. [T]). The straight phase boundaries are 
represented by solid straight lines, but the parabola is distorted into a curved W shape to give 
room between it and the lines for the symbolic labels. Light grey (yellow online) circles label regions, 
darker (blue online) circles label boundaries, and darkest (green online) circles label intersection 
points, x's are extrema of e(u), | and | are u = —1 and 1; one x above the other means the 
minimum and maximum coincide (this occurs on the parabola M = 0). In the three regions in 
which there are minima separated by barriers, a sketch of E vs cos 9 is shown, to define the left, 
middle, and right (L, M, and R) barriers. 

V. PHASE DIAGRAM AT NONZERO TEMPERATURE 

At nonzero temperature, the system will not remain in a well with a very low bar- 
rier. If we assume an Arrhenius-Neel model for the switching rate with a prefactor u, 
an experimental time scale r, and an experimental temperature T, the critical value 
(call it e) of the dimensionless barrier Ae, below which a well is not stable, is given by 
Tuexp(— HqHk MsVe/kgT) = 1; here we use the value 0.03 for e. Only 3 of the regions 
have barriers that can trap the magnetization: | x |x, | x x|, and x| x |. The energy e(u) 
is sketched in Fig. [2] in each of these regions, showing the barriers at the left (L), middle 
(M), and right (R). The zero-T phase boundaries (solid lines) are where L = 0, M = 0, and 
R = 0. The dashed lines where L = e, etc., are the T > phase boundaries. For example, 
in the center of Fig. [2] between the M = and M = e labels, the parallel and precessional 
states are both stable at zero temperature, but at T the precessional state will jump the M 
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barrier and only the parallel state is stable - thus this point is effectively in the P region, 
and the boundary moves down to the dashed curve. 

The actual T > phase diagram is shown in Fig. [3J It is topologically equivalent to the 
cartoon, but because the boundary shifts are very small in places, the topology is easier to 
see in the cartoon. The main result is that the P & PS coexistence region shrinks, and a 
region of telegraph noise appears near its boundary. 




FIG. 3: Nonzero-temperature phase diagram (colored boundaries) - zero-temperature boundaries 
are indicated by solid black lines. Color scale from white to grey (red online) in the PS (precessional 
state) region and the coexisting parallel and precessional (P & PS) region indicates cos 9 for stable 
precession. We have colored part of the P region grey (green online) to indicate that it is a 
physically interesting region - both the M and R barriers (defined in Fig. [2] are lower than the 
critical value e, so the system jumps between the precessional and the parallel state rapidly on the 
experimental time scale - this is telegraph noise. Inset at top shows detail of part of this telegraph 
noise region. Graphical insets show e(u) energy landscapes at the indicated points, as in Fig. [TJ 
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